Derivation of the Ginzburg-Landau equations of a ferromagnetic p-wacve 
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We derive a Ginzburg-Landau free energy for a p-wave ferromagnetic superconductor. The starting 
point is a microscopic Hamiltonian including a spin generalised BCS term and a Heisenberg exchange 
term. We find that coexistence of magnetisation and superconductivity depends on the sign of the 
energy-gradient of the DOS at Fermi level. We also compute the tunneling contribution to the 
Ginzburg-Landau free energy, and find expressions for the spin-currents and Josephson currents 
across a tunneling junction separating two ferromagnetic p-wave superconductors. 
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INTRODUCTION 



In recent years, experimental evidence for materials 
that simultaneously feature superconductivity and ferro- 
magnetism, has transpired from a number of independent 
studiesiiSi^. This has triggered much research into this 
class of materials, widely known as ferromagnetic super- 
conductors (FMSC). There are still many unanswered 
questions concerning the properties of these systems fea- 
turing coexisting spontaneously broken symmetries (bro- 
ken U{\) and 0(3) symmetries). Recently, the supercon- 
ductivity in ZrZn2 has been found to be surface sensitive, 
and might in fact not be a bulk property of the material*. 
A common feature for all the ferromagnetic superconduc- 
tors is that they are superconducting only in the ferro- 
magnetic phase. That is, when pressure is increased in 
such a way as to destroy the ferromagnetic phase transi- 
tion, then superconductivity, if present, is also lost. 

Ferromagnetic superconductivity is of great interest 
from a theoretical point of view. It offers a laboratory for 
studying condensed matter systems with multiple spon- 
taneously broken symmetries, and their interplay. From 
a technological point of view, it is also hoped that het- 
erostructures with FMSC will give rise to novel trans- 
port effects involving both the charge and the spin of 
the electron, which may have the potential for being ex- 
ploited in novel types of devices. Recently, structures 
containing FMSC have been investigated theoretically^. 
This work, however, considers tunneling effects between 
a spin-singlet superonductor in a Fulde-Ferrel-Larkin- 
Ovchinnikov stateS*i coexisting with helimagnetic order. 
On the other hand, little is known of the tunneling ef- 
fects and interplay between ferromagnetism and super- 
conductivity in spin-triplet superconductors coexisting 
with uniform magnetic order. It is the purpose of this 
paper to derive a Ginzburg-Landau model for such a sys- 
tem starting from a reasonable microscopic Hamiltonian, 
and then apply the results to computing tunneling effects 
both in the charge and spin sector, and thus to elucidate 
the interplay between the broken t/(l) and 0(3) symme- 
tries form such systems. Previously, the Bogoliubov-de 



Gennes equations for such systems have been studied in 
some detail^. 

Phenomenological models designed to describe 
FMSC were proposed soon after their experimental 
discoverjiSiifi. These theories can be written down from 
symmetry considerations of the two order parameters in 
the problem, i.e. the theories make no reference to the 
underlying microscopic physics. In ferromagnetic super- 
conductors it is believed that it is the same electrons 
that are responsible for both superconductivity and 
ferromagnetism. In such a scenario, the Cooper-pairs 
must have a magnetic moment, which means that a 
FMSC with uniform magnetic order must be a spin 
triplet superconductor. Microscopic theories attempting 
to explain the positive interaction between electrons 
with aligned spins have traditionally evolved around 
spin mediated interactions, i.e. magnon exchang o^^d^ . 

In Reflisl the phenomenological model of Refs. IgllTol 
was employed in order to find stable solutions and the 
regions were these solutions of the model are stable. The 
phenomenological model used in Refs. I^fiol is given by 



/ (^, M) = fs{^) + fF{M) + fj{^, M) + ]^H\ 



(1) 



were -0 is a three dimensional complex vector describing 
superconducting order, B = {H + AttM) = V x A is the 
magnetic induction, H the external magnetic field ,and 
A is the electromagnetic vector potential. 

The term fsi'ip) describes superconductivity when 
H = M = 0, and is given by 



+ 



(2) 



with 



+ (AV-^ri^-V-^)] + K3{D,iP,y{D,^,), (3) 

where summation over the indices i, j is assumed and the 
symbol 

A = -i^ + 2eA„ (4) 
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of covariant differentiation is introduced. For a discus- 
sion of the superconductivity part of the free energy see 
Refs. ITfillT^ The term /f(M) in ((TJ describes the ferro- 
magnetic ordering of the material and is given by, 

3 , 
fF{M) = CfY,\V,M,\'+af{Tf)M' + -^M\ (5) 
i=i 

This is the standard expression for an isotropic ferro- 
magnet. Furthermore, fj{4>,M) describes the interac- 
tion between the order parameters for ferromagnetism 
and superconductivity, and is given by 

/7(?/',M)=i7oM•(V'xV'*) + (5M2|V'|^ (6) 

In this paper, we will start out with a microscopic 
Hamiltonian that essentially reproduces the Ginzburg- 
Landau model proposed in Ref. ITol As a byproduct of 
such a derivation from a more microscopic model, we ob- 
tain analytic expressions for the coefficients in The 
microscopic theory we start with is not trying to explain 
the mechanism of the non-unitary spin-triplet pairing of 
the electron gas. We simply assume the existence of some 
spin generalised attractive BCS term in the Hamiltonian, 
in addition to a Heisenberg ferromagnetic exchange term. 

As an application of the Ginzburg-Landau equations, 
we look at two FMSC in tunneling contact. We com- 
pute the free energy of the coupling and further find the 
Josephson current and the spin current in the direction 
perpendicular to the plane spanned by the two magneti- 
sation directions. 



II. MODEL 

We start out with a Hamiltonian that is given by three 
terms, one describing free electrons, one BCS term to 



account for spin-triplet superconductivity, and finally a 
Heisenberg ferromagnetic exchange term to account for 
itinerant ferromagnetism. The Hamiltonian of the sys- 
tem is then given by, 



+ 2 X] ''^fe. '='4+9/2, Q'=lfc+q/2,^'=-fe'+<z/2,/3Cfe'+q/2,Q 



where Ck,(7 and ^ annihilates and creates an electron 
in the state (fc, a) respectively. Here, Sk is the usual spin 
operator given by, Sk = J2q,a,i3 '^\+q,a'^aP^q,0 ^^d cr is 
the Pauli matrices. 7(q) = expiq • 5 is the structure 
factor of the underlying lattice, and <5 is a nearest neigh- 
bour vector. A summation over repeated Greek indices 
is implied. 

We are interested in calculating the partition function 
of the system, formally given by 



Z = Tr(e-^^[--^'l) =e-^^, (8) 
were F is the exact free energy of the system. Introducing 
fermion coherent states, £,k,a, and performing a Hubbard- 
Stratonovich decoupling of the two last terms in (0 , we 
arrive at an effective action (in Euclidean time) which 
reads. 



•S'cff 



dr 



a,0 

[^a/3(fe'g)^fc+q/2,/3^-fc+q/2,a + ^M^, q)^lk+q/2,l3^l 



^k+q/2,. 



k.q 



q ^ ' k,k',q 



(9) 



I 

Here, Apa{k, q) and Mq are auxiliary fields in the func- magnetisation, respectively. 

tional integral representation of the partition function. The above action is Gaussian in the fcrmionic fields 

and physically represent the spin-triplet pairing fields and and hence the integral over the fcrmionic fields may be 

performed exactly. We next introduce a Majorana basis 
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now be written, 



, the effective action may 



^''^ = -\ t ^ ^l''^'"^^ - E j^)^"^-" 

- J2 trAt(fe',q)i;;i,A(fe,q)}. (10) 

k,k' ,q 

Here, tr denotes trace over spin indices. Furthermore, 
is a 4 X 4 matrix given by, — Qq^ — S, where 
00 ^ = diag(-iw„ + Sk, -iw„ + Sk, -iw„ - £_fe, -iUn - 
s-k) describes a free electron gas. Here, S describes 
the interaction and pairing correlations, and is explicitly 
given by, 



M -D 
M 



(11) 



where M. = Af-crisa2x2 matrix order parameter 
describing magnetisation and I? = d ■ cr = —iAay is a 
2x2 matrix order parameter describing triplet supercon- 
ductivity. 

The integral over is Gaussian, and hence it may be 
performed analytically. The integral produces a fcrmion 
determinant which may be included in the exponent, i.e. 
in the Ginzburg-Landau free energy, which is given by 



PFGL=-Trlng-^ 



dr 



y — ^MqM_q 



k,k',q 



(12) 



where Tr implies a trace over all variables. The GL free 
energy is defined by 



(13) 



The trace over In^ ^ may formally be rewritten as 
Trlng-^ = Trln^t^^ + Trln(l - ^qS), the first term, 
describing free theory, is neglected in the following. The 
second term is assumed small, and hence we may expand 
the logarithm to obtain 



Trln(l - SoS) ^ - Tr logo's + -Qo'Sgo'E 



-El ~E2~ Ea 



(14) 



The first term in the expansion, Ei = Tr^o^) is zero 
since the Pauli matrices are traceless.The terms E2 and 

axe second order in the ordering fields A^Q(fc, q) and 
Mq, and we now proceed top discuss these terms in turn. 



III. SECOND ORDER TERM 

The second order term in the expansion of the trace is 
given as 



E2 =^ Tr goS^oS = ^ ^ (n| eoS^oS \n) 

n 

= ^ X/ ^Ofei5]fei,fe20ofe25^fe2,fei- (15) 



'»1.'»2 



The last equality comes from inserting completeness rela- 
tions and noting that is local in fc-space. We have in- 
troduced the notation Sofci = (^Sil Qo l^i) and Hki.k2 = 
(fei|S|fe2). After changing variables, } ^ {^+,{2- 
the second order term is 



-E^2 = E { i90,k+q/290,k-q/2 + C.C) Mg ■ M_q 
(90,k+q/29o,k-q/2 + C-c) |dfe,q|^} 



(16) 



where g^ l = (-ia;„ + £fe). 

The trace of the superconducting order parameter is 
given by 



ti:A\k,q)A{k,q) =tr 
=2d 



k.q ' f^k.q- 



(17) 



The complete second order term of the Ginzburg-Landau 
free energy is therefore given by 



^fc,0 _ 1 

Jl{q) P 
1 



-3 {90,k+q/290,k-q/2 + Cc) 



90,k+q/29o,k-q/2 + ^-c) + 2/V \dk,qf } 



(18) 



Now q is assumed to be small, and hence we may expand 
Sk+q/2 « Efe + 9/2 ■ ^ + Oq"^ = Sk + q/2 ■ vp, the 
inverse electron propagator is to first order in q given by 
9o,k+q/2 ~ (-iw„ + efc + q/2 • vf). Hence, we find for 
the part containing magnetisation 



90,k+q/290,k-q/2 + 5fc+q/2.9fe-q/2 



=23? 



1 



{q/2vFr 



{-iuJri + £kf {-iU)n+£k)^^ 



(19) 



when keeping terms to second order in q. In addition, we 
expand the structure factor 7(9) ^ 6 — 2q^, and assume 
for simplicity a cubic lattice in three dimensions. Now, 
we substitute Y^k ~^ ^0 Jl"^ / , where A^o is the 
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density of states at Fermi level and Cc is some cutoff. The part containing derivatives is given by 
The homogeneous part of the magnetisation is given by 



Fl^r. = f E (4^otanh + j^j M,M. 



(20) 



72 



l-tanh^(^) 



— I q^M„M^„ 
36 J J ^ ^ ^ 



(21) 



Similarly, we find for the superconducting order pa- 
rameter, when introducing d^{k,q) = A^i{q)kt'^^ and 
the indices fi and i run from 1 to 3, that the coefficient 
in front of the term involving the superconducting order 
parameter is given by 



q/2 ■ Vp q/2 ■ Vp 



(22) 



Here, we have used vp = vpk. The terms linear in q will 
integrate to zero when integrated over the angles Q. 

We start with the constant term( independent of q) 



An 2Nokikj 



47r iv?-, 



2/37Vo<5,, 



E 



3 ^71+1/2 

n>0 ' 



(23) 



where e is an energy cut-off. In total the constant second 
order term is given by 



2,S 



IE 



2 AN„ 
W ^ ~3~ 



ri>o ' y 



tr ( AA^ 



(24) 

After using the result for the critical temperature^^, the 
term may be written on the form 



2^ 3 Tc 



(25) 



The coefficient for the second order term containing 



derivatives is given by 

a, = -2N,{vp/2)\q,Y.j^^j ,.a^,2 



- 2No{vp/2)'^qiqjl/l5{SijSi„i + SuSj^ + Si„i6ji) 

- ■n/22N(3{vp/2fqiqjl/lb{5ij5hn + SuSj^ + SiraSji) 

,,,3 



^q^qJvlNa|3\{i) 
2407r2 



(26) 



so that the total second order term containing derivatives 
is given by 



7vlNo(3\i3) 
1207r2 

qHrAA^ + {q,A^,) (qjA]^^ + (g^A") 

(27) 



Here we have assumed an isotropic quadratic disper- 
sion relation, hence the coefficients in (O are all the 
same to this approximation. If one introduces an effec- 
tive inverse mass tensor (l/m*)i.j one can use the re- 
sult above by only replacing qi qjvpi j, where vpij = 
l/2[{l/m*)ij + (l/m*)j_i], and remove vp. In this way 
space is no longer isotropic. 



IV. 



THIRD ORDER TERM 



In this section we will find a coupling between the or- 
der parameters for magnetisation and superconductiv- 
ity. The third order term will, however, only be nonzero 
when the system is in a non unitary phase, i.e. when 
d* X d ^ 0. In a ferromagnetic superconductor where 
the same electrons are responsible for superconductivity 
and ferromagnetism, the system, however, has to be in a 
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non-unitary state. If the system is in a unitary state the 
superconducting order parameter is invariant under the 
time inversion operator, and hence there can not be any 
magnetism associated with it. 

The third order term in the expansion of the fermion 
determinant is given by 



3-3 

ki,k2,k3 



Goki Sfci,fc2^0fc2^fe2:fc.?^0 fea^fes^fcl ■ 



(28) 



Via a change of variables, } ^ {I+JS , t} 



{ k+q2/2 



and 



k3 \ 
ki j 



r fc+q3/2 



along with the constraint. 



\k-q2/2 "^^^ fci J \k-q3/2 

Qi + Q2 + Qs = 0, a multipication of the matrices yields 
the following expression for the third order term 

£^3 = 4 ^ {go,k9o,k9a.k + c.c) 



(29) 



were the combination m = id x d* is interpreted as the 
average magnetisation due to the Cooper-pairsi^ii^. The 
prime on the sum in Equation H29|l denotes sumation over 
configurations with the restriction 6{qi + q2 + qs)- 

In the coefficient we have neglected the q dependence, 
since expanding in powers of q amounts to finding deriva- 
tives. The coefficient is given by, when introducing 

53/3! =4K J2 ^ 



<^n,k 



2 I ^2 



lUJn + £k 



AVl ^ ki kj 



=^o/^^./3E;^- (30) 



n>0 



The final summation in the above expression is propor- 
tional to the logarithm of the cut-off frequency. The total 
third order term is thus given by 



"3. 



qi.q2,q3 



3! 



\A -A* 



:Mx. 



(31) 



were 03 — as/p. As the order parameters M and d do 
not couple to second order in the GL free energy, the 
coupling in the third order term is expected to be of 
crucial importance as to whether the two order param- 
eters will coexist in the system or not. Coexistence is 
favoured by the system if there is an energy gain by hav- 
ing M and d finite simultaneously. Hence, coexistence 
of magnetism and superconductivity depends on the sign 
of 03, which again is given by the gradient of the DOS at 
Fermi level. In the system considered here, there will be a 
ferromagnetic coupling between magnetism and the spin 
magnetism of the Cooper-pairs if N^^ > 0, i.e the gradient 



of the DOS should be positive at the Fermi level for the 
spin-up sheet if coexistence of FM and SC is preferred by 
the material. 

Comparing Eq. (|31|l with the model Eq. ©, we find 
the value of the coefficient 70 when we assume that is 
in a p-wave state. The sign dependence of is however 
more general, since a triplet state must necessarily be 
an odd function of k. From H29I) . we observe that in 
an expansion of the DOS around Fermi level, iV(^) ss 
iVo + + CV2A^o + only the odd terms of the 
expansion give a contribution to the coefficient. Hence, 
to lowest order 70 ^ -/Vg. We next determine the fourth 
order terms in the order-paarameter expansion of the GL 
free energy. 



V. FOURTH ORDER TERM 

In this section, we will find three different types terms, 
namely a term involving only magnetisation, one onoly 
involving superconductivity, and finally one term involv- 
ing a coupling between the magnetisation M and the 
superconducting orderparameter d. In total, we will find 
five different independent terms for a p-wave ferromag- 
netic superconductor. The fourth order term of the ex- 
pansion of the fermion determinant is given by 

-£'4 = - E Qoki'^kuk2^0k2^k2,k3 



ki ,fe2 
fc3 , fe4 



GQk3'^k3,k4,GQki'^k4,ki ■ 



(32) 



As before, we change variables and multiply out the ma- 
trices and take the trace. In addition, we make the ap- 
proximation of neglecting the dependence on qi in 50, fci • 



A. Terms containing only magnetic order 
parameter 

The term involving the magnetisation only is given by 



{q.} 



qi) 



^ X! 5o,fc5o,fc.9o,fc5o,fc I • (33) 

Here, the coefficient of the TVfg-factors is given by the 
trace over the electron propagators. 



=5R (5o,fc)' = iVoJR^ / " dd - 



k.UJn 

= ^tanh( — 



1 



1 — tanh 



V ^i^" + ^ 

2 / £m/3 



(34) 
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Thus, the fourth order term involving magnetisation 
only, is given by 

{q^} 

where a^yn = 4!a4m//3. 



The coefficient is given by 



Q4S _ _/Vo_ 

4! ~ 15/3 



^m^^ ^^^^ 



In total, the fourth order term involving superconductiv- 
ity alone is given by. 



B. Terms containing only superconducting order 
parameter 

The term involving the superconducting order param- 
eter only will contain five different terms. Again neglect- 
ing the q dependence in .go.fei we arrive at 



E' 'XS — ^ / ^3 ^\ ^TTl 



E 



Now let the sum over k go over to an integral over ^ and 
the angles 6 and 0. The integral over the angles produces 
Kronecker-(5's in Latin indices. 



(37) 



tvAA 



^ + 2 (trAA^^ ^ 



+ 2tr(^^) 



2tr(^t 



(39) 



where = '^''^^2''^° ■ This part of the free energy thus 
consists of five independent termsi^. 



C. Terms containing magnetic and 
superconducting order parameter 

The fourth order coupling term between the magneti- 
sation and the superconducting order parameter is given 
by 



Sm 



E \ 85o,fc5o,fc5o,/=5o,fc [2 (Mq,d_qJ {Mq,d: 



Qi, 



^Ggo^kgO,k90,kgo^k {dqsd-qj (Mq^M^q^) 

2 {MM {M.Al^) 



^'7iVoC(3)/?3 



(40) 



where the constant in the free energy is a^sm = 
72iVoC(3)/3^/7r^. This coupling term between magnetism 
and superconductivity differs from the fourth order cou- 
pling in Eq. (jB]). In Eq. the coupling term is propor- 
tional to \M ■ ^\ , whereas we find a term proportional 
to {M ■ M) In addition, the coefficient is positive 
rather than having an indefinite sign as commented on 
in Ref|0. We have, however, assumed that the system is 



in a p-wave state. In the a case of a general spin-triplet 
state, both types of fourth order coupling terms between 
magnetism and superconductivity may be present, as also 
seen from Eq. H40I) . 

D. Complete fourth order term 



In total, the fourth order term is thus given by 



7 



V3> 



' + 2 ftr^^V + 2tr ( AAA ( AA^Y + 2tY ( AA^^ ' 



F4 = 1/4! J2 

{Qi} 

+a4s - tr AA^ 
-2tr (^AA^^ (A4t) 

All in all the Ginzburg-Landau free energy for a ferromagnetic superconductor is given by 
Fgl = y" d^r I tr AA^ + ^ (d^ tr AA^ + D^Af^.DjA]^ + D,A^,^D^AJ, 



(41) 



"45 

4! 

Q:4m 



2(tr^^)- tr^A^ +2tr(^^)(^A^)* + 2tr(^^)2-2tr(^^^)(^^)* 



(M • M) 



2 Q4m 
4! 



(42) 



where D = V + 2ieA, D = V + ieA and A is the elec- 
tromagnetic vector potential. 

VI. TUNNELING 

As an application of the Ginzburg-Landau theory de- 
rived above, we will consider tunneling between two fer- 
romagnetic p-wave superconductors. The Hamiltonian 
we use is the common choice when studying tunneling 
between systems in equilibrium 

H = Hl + Hr + Ht, (43) 

where Hj^^ji^ is given by iQ and 

= E + T*^'dlpc,,»] ■ (44) 

h.p 

For a general Hamiltonian, it is straightforward to 
show that the part of the Ginzburg-Landau free energy 
that contains the tunneling elements, is given by 

Fj = -]-TT\n{l-gi^Tg^T), (45) 
P 

where ^l(r) are the Green's functions for the left(right) 
subsystem and 7" is a tunneling matrix, 



with 







-fc,— p 



k.p h.p 

r-pi\ rpil 

k,p k,p 



(46) 



(47) 



Furthermore, we will assume for simplicity that the 
magnetisation of our system is weak, hence wc assume 
that the magnetisation is homogeneous. Introducing 



and choosing the 



the basis cpl = ?I,|,?I,x,?-fc,T,C-fca 
quantisation axis along the magnetisation vector at each 
side of the junction, we have the following effective ac- 
tion, 



ScS = S + S ' + S 
and the inverse Green's functions are given by 




(48) 



(49) 



direc- 



tion of the magnetisation does not enter in the Greens 
functions for the isolated systems on the left and right, 
the angle between the magnetisations, and hence the 
quantisation axes, only enters through the tunneling el- 
ements since we have applied different rotations on the 
left and right. This difference shows in the part of the 
Hamiltonian where the left and right systems couple. 

We now want to calculate the Ginzburg-Landau free 
energy for the junction, Fj. To this end, we expand l|45|l 
to the lowest order in the tunneling elements. 



1 



Fj^-Tv[gLTgRT] 



(50) 



Next, we use the Dyson equation to find an expansion of 
g i.e. g fa go + g^'Sgo + ... were 



g 

-g* 



(51) 
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and 



S = 



A 

At 



(52) 



Since we are assuming that the magnetisation is homo- 
geneous, Qo is local in fc— space, S however is nonlocal 
as we do not want to impose any restriction on the su- 
perconducting pairing state. 

From the free energy of the junction we can find the 
current of various quantities associated with tunneling. 
To show this we introduce a Hamiltonian H {q,p) which 
is a function of the canonically conjugate variables q,p^ 
and recall that the velocity operator is given by g = ^ . 

Consider now the quantity || = -i^ In [Tre"'^^] = 



J-Tr 



dH^-0H 



(q). Hence we are easily able to to 

find the current over the junction of for instance elec- 
trons. Since [0, N] — i, the current is than given by 

I J = — e (^n'^ = e-^Fj and the contribution to the spin 
current coming from tunneling across the barrier, is given 
by (^Sn^ = —^^b-§qFj, where [9, Sn] = i, and 9 is the an- 
gle in the plane perpendicular to the direction fi. 

A. Ferromagnetic case 

We start with the term corresponding to a junction 
between two ferromagnetic metals, 

^J^^ ""^ {^okTk,pQQpTl p 



i J2 tr{gfeTfc,pgpTt p 



(53) 



We observe that upon letting k ^ —k and p p the 
second term above is the complex conjugate of the first 
one, hence the free energy is two times the real part of 
the first term. 



F 



(0) 



IS" 



9k9p 



'9k 9p 



rp+ + 

^k,p 



9k 9p 



T, 



k.p 





2 






k,p 


+ 9 k 9p 


k,p 





(54) 



Here -I- and - means parallel or anti parallel to the mag- 
netisation respectively. The angle between the magneti- 
sations on the left and right enters through the tunneling 



elements 



T, 



a/3 



= 1/2 [1 + af3M^ ■ M'^j \Tk,p\ , when 

we assume that the spin in preserved across the tunnel 
barrier. The potential difference between the left and 
right subsystems are taken care of through the tunneling 
elements i.e. we assume that Tk,p is nonlocal in w-space 
with a boson frequency lo^, later we do an analytical 



continuation Ui, 
functions is, 



9k9^p = 



eV -\- \5. The structure of the Green's 



(55) 



-i(w„ 



-lUJr, 



were £ka = £k — aM^ and similarly for the right subsys- 
tem. The zero'th order term of the free energy may be 
written in two parts, one part independent of the direc- 
tion of the magnetisations and one part proportional to 
the dot product of the directions of the magnetisations, 
that is, to the cosine of the angle between the magneti- 



sations, F 



(0) 



ML . MRpf]^. With 



k,p 

2 



^ [5fc ffp + 9k 9p + 5fc .9p + 9k 9't\ \Tk.:p\ 



d& / d^pA^(a)^(^p) 



/(& ^ 


~ M^) - 


/(Cp- 








f - 


-M^ 


/(a 


^ M^) - 


- fi^P - 


- M«) 




- £,k - 




-M^ 


/(a 


- M^) - 


- fi^P - 


l-M^) 


eV ^ 


-ik- ip 




- M« 




+ M^) - 


- fiip ^ 


- M^) 


eV- 


^^k-^p 







(56) 



and 



-AO) 

J,M 



1 



2^ [9k 9p + 9k 9p - 9k 9p - 9k ff; 



ITI- 





- M^) ~ 


fi^p- 




eV + 


^k - £,p' 


1- - 


-M^ 


f{ik 


-M^) - 


- /(Cp - 




eV ^ 


-£.k- £,p 




-M« 


fi^k 


- M^) - 


- /(Cp - 


V M^) 


eV ^ 


-ik-s.p 




-M« 


fi^k 


+ M^) - 


- /(Cp - 


-M«) 


eV- 


\^ £,k - £,p 







(57) 



Here N{S) is the density of states and f{x) = 1/(1-1- 
e'^^) is the Fermi-Dirac distribution function. The spin 

current (^Sn^, were n denotes the direction perpendicular 

to the plane spanned by the magnetisation vectors, is now 
easily found by taking the derivative with respect to the 
angle between the magnetisations, i.e. 



Sr 



--f^sF 



(0) 



sm( 



(58) 



This is precisely the same result as found previously in 
Ref. 113 via a different route. 
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B. Ferromagnetic superconducting case 

In this subsection, we calculate the free energy which 
is first order in the superconducting gap function on both 
the left and right side, i.e. we consider the ferromagnetic 
superconducting state. We want to find the Josephson 

current and also the two particle contribution to (^Sn^- 

The free energy which can give rise to Josephson current 
is, 

fP -1 Tr [goL'SLQoLTgoRi:RgoRT] 

(59) 

We now change variables, t ^ and K P+Vji 

^ ' k' k — q/2 p' p— g/2 

were q{q) is the centre of mass momentum of the Cooper- 



pairs on the left (right), and hence it is small. Expanding 
in q{q) amounts to an finding derivatives, we are not in- 
terested in derivatives and neglect the q{q) dependence 
in all terms except S, since we want to have the possi- 
bility ^ 0. With these assumptions the free energy 

reads, 

=^Etr [g;.A(fc,q)g*_^T*_^,_^g*pAt(p,-g)Tt p 
+g-fe A(fe, -q)gfcTfc,pgpA(p, q)g;T*_fc _p . 

(60) 

Here we observe that upon changing the sign of all im- 
pulses in the second term and using that A(— fc) = 
— A(fc) for a triplet superconductor, the second term is 
just the complex conjugate of the first one. Hence, the 
free energy is just two times the real part of the first 
term. After some straightforward algebra, we obtain 



fT = -\t\'Y. 



{( 



1 + • M 



+ ^ii^^// + A|i(F;/ + ) 1 cos A0 



1 -l-i TT'-ii 



A}|;i}| + A||;i|| + A}|(;i}| + /i||))sinA0 



1 - • M 



J\_ ^ ^ J_ "J" "I" 



"^TT TT ^ ii ii Ti*> it 



cos A(/) 
sin A0 



a{|(f;/ + F/; ) - a||(f;; + f;; ) + ^-(f// + f//) - A|i(F;/ + f- ) j cos a^ 

+ (A\\{h\\ + h\\) ^{{(/ill + h^^ + A\\{h\\ + h\\) A\\{h\\ + h^S) ^ii^H } 
were we have introduced the notation 

/n„>o 47r 



Jnk>o Jo,, 



R I 



(62) 



F 



A7 



Q/3 



dCfed^, 



(A^(-a + c^M^) + N{^k + PM^)) {N{^p + AAf«) + N{-^p + 7M«)) /(^fc) - /(^p) 



and 



4 Uk - ^A/M ( Cp + 



eV^ + a - Cp 



(63) 



^ T / da [/(a) - /(a 



N{^k + a M^ + eF) {Nj^k + iM^) + jV(-a + \M^)) 
ik + + ev) (a - V^^"" 

N{ik + /3M^) (iV(a + AAf ^ + e\^) + jV(-a + -fM^ - eV)) 

(ik ~ ^M^) (Cfc + V*^'' + 



(64) 



In addition, we have assumed that — (j)^^ — At/) is independent of the spin indices, this can be done since any 
coupling between the different components of the superconducting order parameter will give a phase lockinglS. Further 
since Tf, _p = OiS the sum over k and p i Equation H59|) may be restricted to positive values only, and hence there 
are no problems involved in splitting the superconducting orderparameter into an amplitude and a phase in Equation 
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(|61|l . The Josephson current is now found by simply taking the derivative of the free energy with respect to the phase 
difference A(p. This leads to the floowing expression for the Josephson current 



e\T\ 



1 + ■ M 



")[( 



/iTTpTT 



n 



sin A4> 



1 - ■ M 



j_ 



/iTTpTT _ A^HF 



sin A(/) 



- ( A\\h\\ + A\\h\\ - A\[{h\\ + /iji)) cos A0 



TT TT 



pTT^ 



1 ii 



Similarly, we find the following expression for the two-particle contribution to the spin current 



cosA0 I . (65) 



Sn)^ = -\T\ 



X{sin.[( 



/iTT pTT 

21. '|- J. 'j" 



"^ii^ii 



iTi/pTi 
Hi^^Ti 



it 
iT 



it 
Ti 



iT^ 
Ti > 



— /liipii _ 4TTpTT 
^TT-^TT "^ii^ii 



cos A0 



(^Ti^^iY + ^tY-* ^ ^}|(-F|Y 



cos 6* 



"ti 
i^lY ) 



ijiTN _ ^iiLii _ /iTTlTT 
"'Ti'' "^ii^ii 



sin A(/) 



piT^ 
■^ii ' 



\ ii / pii 
^i^-^Ti 



pii ■ 
-^iT • 



cos A0 



+ (^Ti(/YT + ^U) - ^TT(^n + ^Tt) + ^ii(^U + Mi) - ^Ti(Mi + M})) H } 



(66) 



The spin current is not well defined in the limit when 
— > 0, since we are calculating the current along the di- 
rection perpendicular to the plane spanned by the two 
magnetisation directions. We observe that the total spin 
current is the sum of Eq. H58|) and (|66|l . Hence, Eq. 
(|66|l is a contribution originating with the interplay be- 
tween superconductivity and magnetism. It is seen to 
disappear when superconductivity is lost. Moreover, we 
observe that there is a term in equation H66() that is pro- 
portional to COS0. In Josephson currents there also exist 
a term that is proportional to cos A0, which however is 
dissipative and vanishes if the voltage across the junction 
is set to zero. In the case of the spin current, we find that 
the condition that must be fulfilled for the cosine term 
to go to zero, is given by A^^ = A^^ , — 

and eV — 0. Furthermore, it is seen that the cosine term 
is associated with flipping one electron from (spin)state 
a on the left to —a on the right, hence the cosine term 
is proportional to | A^^ | | A^^^ | and similarly with the 
superscripts interchanged. 



between the two order parameters is a third order term, 
this term is only non zero when the material is in a non 
unitary superconducting state i.e. d* x d ^ Q. From 
an exact calculation of the coefficient of the third order 
term we find that coexistence of ferromagnetism and su- 
perconductivity is enhanced if the gradient of the density 



of states, ^N{C) 



is positive at Fermi level. 



We have also computed the tunneling contribution to 
the Ginzburg-Landau free energy. From this, we have 
found expressions for spin- and charge currentws in the 
spin- and charge channels across a tunneling junction sep- 
arating two ferromagnetic p-wave superconductors. 

Superconductivity coexisting with ferromagnetism can 
be triggered by the magnetisation as a result of the pres- 
ence of a third order term^^^, i.e. superconductivity ap- 
pears at a higher temperature than if magnetisation were 
not present. The sign of the third order term is given by 
the the sign of the gradient of the DOS at Fermi level. 
Since the Fermi level may be tuned by applying pressure, 
it may be possible to change the sign of the third order 
term by applying pressure. 



VII. 



CONCLUDING REMARKS 



We have derived a Ginzburg-Landau functional from a 
microscopic Hamiltonian consisting of three terms, free 
fermigas, an Heisenberg term and a spin generalised BCS 
term. We find two order parameters, local magnetisation 
and the superconducting gap. We expand to fourth or- 
der in the order parameters. The lowest order coupling 
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